Astronomy & Astrophysics manuscript no. aal7896-ll 


©ESO 2012 


February 14, 2012 





Rubidium, zirconium, and litliium production in intermediate-mass 

asymptotic giant brancli stars 

M. A. van Raai\ M. Lugaro^, A.I. Karakas^, D.A. Garcia-Hemandez'*'^ , and D. Yong^ 

' Sterrekundig Instituut, University of Utrecht, Postbus 80000 3508 TA Utrecht, The Netherlands, e-mail: markvanraaiOgmail . com 
^ Monash Centre for Astrophysics (MoCA), School of Mathematical Sciences, Monash University, Clayton Victoria 3800, Australia, 

e-mail: maria . lugaro@monash . edu 
^ Research School of Astronomy and Astrophysics, Mt. Stromlo Observatory, Cotter Rd., Weston, ACT 2611 Australia, e-mail: 
akarakasOmso . anu . edu . au , yong@mso . anu . edu . au 
' Instituto de Astrofisica de Canarias, C/ Via Lactea s/n, 38200 La Laguna (Tenerife), Spain, e-mail: agarciaOiac . es 

D , ^ Departamento de Astrofisica, Universidad de La Laguna (ULL), E-38205 La Laguna, Spain 

[-■h ■ 

_ . Received / Accepted 

m. 

ABSTRACT 



(N 

o 



Context. A recent survey of a large sample of Galactic intermediate-mass (>3 Mq) asymptotic giant branch (AGB) stars shows 
that they exhibit large overabundances of rubidium (Rb) up to 100-1000 times solar In contrast, zirconium (Zr) is not enriched 
compared to the solar abundances in these stars. These observations set constraints on our theoretical notion of the slow neutron 
(— I _ capture process (i- process) that occurs inside intermediate-mass AGB stars. Lithium (Li) abundances are also reported for these stars. 

^ ^ In intermediate-mass AGB stars, Li can be produced by proton captures occuring at the base of the convective envelope. For this 

I ■ reason the observations of Rb, Zr, and Li set complementary constraints on different processes occurring in the same stars. 

O ' /A/'ms. We present predictions for the abundances of Rb, Zr, and Li as computed for the first time simultaneously in intermediate-mass 

AGB star models and compare them to the current observational constraints. 
C/5 ■ Methods. We calculate the Rb, Zr, and Li surface abundances for stellar models with masses between 3 and 6.5 Mq and metallicities 

■ between 0.02 and 0.004. 

Results. We find that the Rb abundance increases with increasing stellar mass, as is inferred from observations but we are unable to 
match the highest observed [Rb/Fe] abundances. Variations of the reaction rates of the neutron-capture cross sections involved with 
Rb production and the rate of the '^Ne(a, «)^^Mg reaction, responsible for neutron production inside these stars, yields only modest 
variations in the surface Rb content of ~ 0.3 dex. Inclusion of a partial mixing zone (PMZ) to activate the "C(q-, n)'*0 reaction as 
an additional neutron source yields significant enhancements in the Rb abundance. However this leads to Zr abundances that exceed 
the upper limits of the current observational constraints. If the third dredge-up (TDU) efficiency remains as high during the final 
stages of AGB evolution as during the earlier stages, we can match the lowest values of the observed Rb abundance range. We predict 
large variations in the Li abundance, which are observed. Finally, the predicted Rb production increases with decreasing metallicity, 
fS) ■ in qualitative agreement with observations of Magellanic Cloud AGB stars. However stellar models of Z = 0.008 and Z = 0.004 

I intermediate-mass AGB stars do not produce enough Rb to match the observed abundances. 
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1 . Introduction Zirconium (Zr) on the other hand is a typical element made by 

^ ■ the s process, with ~ 15-30% of its solar abundance ascribed to 

[ The elements heavier than iron are produced almost entirely by the r process. The r process is believed to occur in supernovae 

neutron-capture nucleosynthesis via the i(low) and the r(apid) and/or neutron star mergers (lMeveilll994l) . while the s-process 

processes, which are each responsible for roughly half of the elements from Sr to Pb are produced in asymptotic giant branch 

cosmic abundances of these elements (Kappeler et al. 1989). (AGB) stars (Gallino et al. 1998). 

During the s process timescales for neutron captures on unstable The AGB is the final nuclear burning stage in the evolution 
isotopes are typically lower than the decay timescales, which im- of stars with initial masses in the range 0.8 to 8 MgO. During 
plies neutron densities A^„ ~ 10^ cm I Conversely, during the r the AGB the H-burning shell is responsible for the nuclear en- 
process timescales for neutron captures on unstable isotopes are ergy production most of the time, but every 10^"^ years or so the 
higher than decay timescales, which implies N„ > 10^" cm^^ He shell ignites, driving convection throughout most of the He 
The element rubidium (Rb) is an example of an element that intershell, the He-rich region located between the H- and the He- 
may have a significant r-process component, where -20-80% burning shells. This is known as a thermal pulse (TP). After a TP 
of its solar abundance has been ascribed to the r process. The has occuiTed, at which time the H-burning shell lies dormant, the 
exact value of this component depends on the i-process model convective envelope extends inward and may reach the He inter- 
used to calculate it, which can be parametric ( Arlandini et al. shell. This is known as the third dredge-up (TDU) and can oc- 
Il999t ISimmerer et"an i2004). stellar ( Arlandini et al. il 19991). or a 

chemical evolution model of the Galaxy (TravagUo et al.||2004- ' Hereafter we refer to AGB stars with initial masses < 3 Mq as 

The r process contribution also depends on the uncertainties i n "low-mass" AGB stars, and to AGB stars with initial masses > 3 Mq as 

the solar abundances and nuclear physics inputs (lGorielvlll999h . "intermediate-mass" AGB stars. 
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cur periodically after each TP. The TDU mixes nuclear burning 
products from the He and H-burning shells to the stellar surface, 
increasing the abundance of carbon and in some cases leading to 
C>0 at the stellar surface ( Herwigll2005b iKarakas & Lattanziol 
l2007HLugaro & Chieffilllom 

In the He intershell of AGB stars, neutrons for the s pro- 
cess can be released by the ^^C(a, n)'*0 reaction during the 
periods in-between TPs (interpulse) and by the ^^Ne(a, n)^^Mg 
reaction in the convective TPs. The '^C neutron source oper- 
ates from T > 0.9 X 10^ K and provides N„ ~ 10^ cm"^ 
In order to produce enough '^C to match the enhancements of 
the s-process elements observed in AGB stars, it is assumed 
that some protons are mixed down into the He intershell at the 
deepest inward extent at each TDU epi sode (Gallino et al. 199^ 
iGorielv & Mo wlavi"2000' 'Busso et al 2001). We refer to the re- 
gion of the He-intershell that has undergone some partial mixing 
of protons as the "partially mixed zone" (PMZ). The protons 
are captured by the abundant '^C producing '^C and '"^N, which 
allows f or the ^^C(g, »')'^0 reaction to be efficie ntl y activated 
(Hollow eiretanil990l: iGorielv & Mowlavill2000l: Eugaro et alJ 
12003). The ^^Ne neutron source requires T > 3.5 x 10^ K to be 
significantly activated and produces A^„ « 10'° - 10" cm"-' in 
low-mass AGB stars and up to A^„ x lO'-' cm"-' in intermediate- 
mass AGB stars (see Sect. [Hi. This is because temperatures in 
excess of 3 x 10^ K are more e asily reach ed during the TPs 
in intermediate-mass AGB stars (Ilbenlll9'75h . The '^C neutron 
source is active for a much longer period of time (~ 10^ yr) than 
the ^^Ne neutron source (~ 1 yr) and it produces a higher total 
amount of neutrons, even though with a lower neutron density. 

The amount of Rb produced during the s process depends on 
the probability of the two unstable nuclei **^Kr and ^*Rb captur- 
ing a neutron before decaying and acting as "branching points". 
The probabil ity of this happen ing depends on the local neutron 
density (,Beer & Macklinlll989l) . Figure [T] shows a section of the 
chart of the nuclides from Kr to Zr. iBeej (Il991h reported that, 
independently of the temperature, when ^"^Kr captures a neutron, 
50% of the flux proceeds to the ground state of **^Kr (with half- 
life of 3934.4 days) and the other 50% proceeds to the metastable 
state of ^^Kr (with half-life of 4.480 hours). Of the metastable 
^^Kr, 20% decays to its ground state while the remaining 80% 
decays into ^^Rb. Effectively, 40% of ^*Ki+n results in the pro- 
duction of "^^Rb, whereas the other 60% results in the production 
of the ^'^Kr ground state. The long decay time of the **^Kr ground 
state allows it to capture another neutron to produce **^Kr instead 
of **^Rb, provided the neutron density is higher than a 5 x 10^ 



cm" . Kj can also capture a neutron, producing Kr which 
quickly decays (with a half-life of 76.3 minutes) into '^^Rb. The 
other unstable isotope in Fig.[T]is ''^Rb (with a half-life of 18.63 
days). This branching point produces ^^Rb directly, with » 50% 
of the fl ux going to ^^Rb if the n eutron density is > 10'° cm"^ 
(see also lLugaro & Chieffill20TTI) . 

In summary, if the branching points at **^Kr and ^''Rb are 
open ^^Rb is produced. This nucleus has a magic number of 
neutrons of 50 and a low Maxwellian-averaged neutron cap- 
ture cross-section of 15.7 mbarn at a thermal ene rgy of 30 keV, 
as compared to 234 for ^^Rb dHeil et al.1 [20081) . This means 
that if ^^Rb is produced, it tends to accumulate. Therefore, the 
^^Rb/^^Rb isotopic ratios is a direct probe of the neutron density 
at the production site. It is not possible to distinguish i ndivid- 
ual '^^Rb and '^^Rb from stellar spectr ( Garcfa-Hernan dez et alj 
120061: iLambert & Lucklll976l) . However, it is possible to deter- 
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Fig. 1. The .s-process path in the region of interest in the chart 
of the nuclides. Isotopes represented by white boxes are stable, 
isotopes represented by black boxes are unstable. The black ar- 
rows indicate the ^-process path. The solid arrows indicate the 
.s-process path branches specific to low neutron densities. The 
dashed arrows indicate the .s-process path branches followed un- 
der high neutron densities. The rectangular box encircles iso- 
topes with magic neutron number = 50. These isotopes have 
relatively small neutron capture cross sections. 



mine the elemental Rb abundance in relation to other neighbour- 
ing 5-process elements, for example Sr, Y, or Zr This ratio is 
usually expressed relative to solar and in a logarithmic scale: 
[Rb/Zr] or [Rb/Srfl A positive/negative [Rb/Zr] or [Rb/Sr] ratio 
implies a high/low neutron density at the ^-process production 
site. This fact has already been used to conclude that the neu- 
tron density in M, MS, and S stars must be low, and that the 
'-'C(a , «)^^0 reaction mus t be the main neutron source in these 
stars dLambert et alJ[T995l) . Theoretical models predict that the 
activation of the ^^Ne neutron source and the maximum value of 
the neutron density are expected to be related to the AGB stel- 
lar mass. For this reason the [Rb/Zr] ratio in AGB stars is also 
predicted to b e positively correl ated with the initial mass. Using 
this argument lAbia et alJ (1200 Ih concluded that C stars must be 
AGB stars of relatively low mass, < 3 Mq. 

Intermediate-mass AGB stars over > 4 Mq also experience 
"hot bottom burning" (HBB). This occurs during the quiescent 
interpulse phase when the convective envelope extends so far 
inward that it touches the H-burning shell and allows proton- 
capture reactions to take place at the bottom of the envelope. The 
CN cycle is activated, which converts '^C into '"^N and keeping 
C<0 in the AGB envelope by reducing the C content. Another 
consequence of HBB is the production of lithium (Li) i n form 
of ^Li via the Cameron-Fowler mechanism (Cameron & TruranI 
ll977l : lBoothrovd et alJll993l) . This happens via activation of the 
'*He(^He,7)^Be reaction at the base of the envelope, followed 
by transport of ^Be towards the stellar surface. There, ^Be turns 
into ^Li via electron captures while the destruction of ^Li by 
proton capture is inhibited by the low temperatures. Models of 
intermediate-mass AGB stars find that the stars become Li-rich 
during the first ~ 10 TPs (Sect. 14.4b . Hence, an enhanced Li 



^ Though they ca n be measured in the interstellar medium 
dPederman et alj2004) . 



[x/Y] = iog,o(x/r) - iogio(Xo/y0). 
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abundance c an be interpreted as a sig nature of the AGB star be- 
ing massive (ISmith & Lamberllll9 90rl. 

Studies of intermediate-mass AGB stars in our Galaxy and 
in the Mage l lanic Clouds (OH/IR stars, iGarcia-Hernandez et alJ 
120061 120071 [2009) have yielded the first abundances of Rb, 
Zr, and Li in these stars. The aims of this paper are twofold. 
First we compare our predictions of Rb and Zr production in 
intermediate-mass AGB stars to the observations in order to re- 
confirm that the ^^Ne(a,n)^^Mg reaction is the main neutron 
source. Second, we try to use the abundances of Rb, Zr, and Li 
and their observational trends (e.g., Rb increases with increas- 
ing stellar mass and/or decreasing metallicity) to constrain the 
uncertainties in the stellar modelling of intermediate-mass AGB 
stars. From our detailed comparison we will attempt to define 
which stellar modelling uncertainties (or other uncertainties e.g., 
model atmospheres) play the largest role in the apparent mis- 
match between observations and theory. We will also discuss 
possible directions for future work to study these problems in 
more detail. 



2. Observations 

The study of the brightest OH/IR AGB stars in our Galaxy and in 
the Magellanic Clouds by ' Garcia-Hernandez et a (120061120071 
[2009) showed that these stars are extremely Rb rich but not Zr 
rich. This suggests that the ^^Ne(a, «)^^Mg reaction must be the 
dominant neutron source in intermediate-mass AGB stars. 

The mass of the observed stars can be estimated from 
the observed Mira pulsation period if the distance to the star 
is known (see, e.g.. Wood et al. 19_83). For the stars in the 
Magellan ic Clouds this comparison indi cated masses of at least 
6-7 Mq jGarcia-Herna ndez et al.l l2 009l). For the observations 
of Galactic stars of iGaicfa-Hernan dez et alj (l2006i l2007h the 
distance is unknown. However, the Mira pulsation period for 
these stars shows a correlation with the OH expansion veloc- 
ity (vexp(OH)). Independently for this sample, the correspon- 
dence of the Vexp(OH) with mass is obtained using the loca- 
tion of the observed stars with respect to the Galactic plane . 
( See Sect. 5.2 and Sect. 5.3 of iGarcia-Hernandez etal] 120071 
for a full discussion on the estimation of mass from Vexp(OH)). 
From the relationship b e tween Vexp(OH) and the [Rb/Fe] ratio 
iGarcia-Hernandez et al.l ( l2006h estimated that the [Rb/Fe] ra- 
tio increases with the initi al stellar mass (see their Fig. 2 ). For 
the same sample of stars ' Garcia-Hernandez et al.l (2007) also 
measured the Li abundances and found a large variations with 
log(Li/H) + 12 ranging from -1.0 to 2.6. 

In Table[T]we summarize the Rb and Zr abundances observed 
in low-mass and intermediate-mass Galactic and Magellanic 
Cloud AGB stars. The putative low-mass stars in the Magellanic 
Cloud reported in Table [U indicate d in brackets are for the 
stars observed by iPlez et al" I (fT99l in the Small Magellanic 
Cloud and for the stars classified as "HBB-AGB" in Table 1 of 
[Garcia-Hernandez et al. (2009). These stars are poor in R b, but 
rich in Li, with log(Li/H) + 12 between -1-1.9 and -1-3.5 (jPlez et alj 
Il993t) . This composition indicates the activation of HBB but 
not of the ^^Ne neutron source. These stars may be more ap- 
propriately considered as a link between low-mass Rb-poor and 
intermediate-mass Rb-rich AGB stars in the Magellanic Clouds. 
We will discuss these stars a gainst mo del predictions in Sect. 14.41 
Table 1 of iGarcfa-Hernandez et al.l (12009i) also reported two 



"Non-HBB-AGB" stars, which have abundances consistent with 
the values reported in Table [T] for putative low-mass AGB in 
the Magellanic Clouds, except that one has [Rb/Zr]>0.3. We 
note that the distance to the Magellanic Clouds makes even 
these brightest of AGB stars faint, thus they are very challeng- 
ing to observe with even the VLT. This means that there are 
only a few observations o f intermediate-mass AGB stars in the 
Magellanic Clouds (e. g.. IWood et al.l 119831: [Smith & Lamberll 
1990; Garcia-Hernande z et al j l2009h . Note that in the sample of 
Garcia-Hernandez et al. (2009*) there are Rb abundances for only 
four stars in the Large Magellanic Cloud and one in the Small 
Magellanic Cloud. These few stars may represent the first truly 
massive extragalactic AGB stars, as indicated by their strong Rb 
I lines. On the other hand, the Galactic sample of intermediate- 
mass AGB stars is much more extensive (120 sources were 



studied bv lGarc fa-Hernan dez et al.l[2006l[2007[) . The main trend 
derived by [Garcia-Hernan dez et al.l (I2009h by comparing the 



* Also -10% of C-rich AGB stars are Li rich dAbia etalJll"993h . 
which points to the occurrence of some kind of "extra mixing" below 
the formal border of the convective envelope in low-mass AGB stars. 



Galactic to the Magellanic Cloud stars is that higher Rb produc- 
tion is observed in lower metallicity environments, with roughly 
[Rb/Fe] >3 in the Large Magellanic Cloud, and [Rb/Fe] <3 in 
their Galactic counterparts. 

Overall, in low-mass stars the [Rb/Zr] abundance is nega- 
tive, while in intermediate-mass stars it is positive. This is in 
qualitative agreement with the current theoretical scenario for 
the s process in AGB stars. Note that the [Rb/Fe] observed 
in intermediate-mass AGB stars has a large uncertainty of + 
0.8 dex, which reflects the sensitivity of the Rb abundance 
to changes in the atmosph eric parameters (Tptf, gravity , etc.) 
adopted for the modelling (IGarcia-Hernandez et al.[ |2006). The 
observational uncertainties allow for an upper limit of [Zr/Fe] 
>0.5 for Galactic AGB stars, and [Zr/Fe] >0.3 for the Magellanic 
Cloud intermediate-mass AGB stars. By error propagation, the 
[Rb/Zr] ratios have maximum error bars of +1.0 dex. We refer 
the reader to the literature sources indicated in Table [T] for more 
information on the abundance derivations and their related error 
bars. 

Serious problems are present in current models of the atmo- 
spheres of AGB stars. These models are performed in ID and 
do not include important effects such as the presence of a cir- 
cumstellar dust envelope and dust formation. Furthermore, abun- 
dance corrections resulting from non local thermodynamic equi- 
librium (NLTE) need to be taken into account. To our knowl- 
edge, such abundance corrections for Rb do not exist in the lit- 
erature, nevertheless we offer a few comments on NLTE. As for 
most, if not all, NLTE calculations, missing atomic data may be 
an issue for Rb. Although NLTE Rb corrections have not yet 
been established, some expected behaviors may be anticipated 
based on N LTE correc tions for other alkali metals such as Li, Na, 
and K (e.g. lPlez et an[T9 93: Asplund 2005; Lind et al. 201 1). In 
fact, for Li and Na lLind et al.l ( 1201 II) noted "many striking simi- 
larities between the two elements, especially in the shape of the 
abundance correction curves." For Na, the maximum NLTE cor- 
rection is about 0.70 dex. If the NLTE Rb corrections mimic 
the behavior for Li and Na, we may therefore expect (a) res- 
onant scattering to be the dominant NLTE mechanism, (b) the 
NLTE correction to be negative, and (c) the magnitude of the 
NLTE correction to reach a maximum for saturated lines. The 
observed Rb abundances in intermediate-mass AGB stars come 
from strongly saturated Rb I lines and therefore, the NLTE Rb 
abundances would be expected to be smaller, although the exact 
magnitude of such NLTE Rb correct ions have yet to be estab- 
Ushed. We note that lPlez et all d 19931) concluded that NLTE Rb 
effects in the stars they analysed are likely to be small. 
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With regard to the 3D hydrodynamical models, the 
temperature-sensitive Rb I resonance Hnes originate not only 
from the photosphere of these stars but also in the outer nonstatic 
layers of the stellar atmosphere and in the expanding circumstel- 
lar shell (Garcfa- Hernandez et alJ '2006). Recent years have seen 
a rapid development of numerical 3D radiative hy drodynamic 
simulations of stellar surface convection (see e.g., ICoUet et"aD 
I2OII . and references therein). In these simulations, gas flows in 
the highly stratified outer layers of stars are modelled by solv- 
ing the hydrodynamic equations of mass, momentum, and en- 
ergy conservation and accounting for the energy exchanges be- 
tween matter and radiation via radiative transfer In this frame- 
work, convective motions arise and self-organize naturally with- 
out the need to introduce adjustable parameters. One of the main 
goals of 3D modelling of stellar surface convection is therefore 
to provide a more realistic description of the physical structure 
of late-type stellar atmospheres. From the point of view of spec- 
troscopic analyses, the structural differences between 3D and 
hydrostatic (ID) models can lead to large differences in terms 
of derived abundances whenever temperature-sensitive spectral 
features are used tAsplund.2005i) . In short, 3D modelling shows 
tiiat 3D eff'ects are more important at low metallicity and low 
gravities. This means that the Rb abundances obtained from the 
sensitive-temperature Rb I resonance lines would be higher in 
the hydrostatic models compared with the 3D hydrodynamical 
ones, as a consequence of the higher temperature in the static 
case, since the strength of t he Rb I resonance lines incre ases 
with decreasing temperature (iGarcfa-H ernandez et al.' '2006"). We 
would expect drastic 3D effects in intermediate-mass AGB stars 
at low metallicity. In summary, possible NLTE and 3D Rb cor- 
rections could decrease the Rb abundances derived from hydro- 
static LTE models by (at least) one order of magnitude. However, 
these speculations need to be confirmed by detailed 3D hydro- 
dynamical simulations and NLTE calculations. 



3. Stellar models 

We ca lculate the nucle osynthesis with a detailed post-processing 
code ('Cannon] Il993h for which the stellar evolution inputs 
were calculated beforehand . The stellar evolution models are 
from iKarakas & Lattanzi 3 (I2007h and were calculated with 
the Monash version of the Mount S tromlo Stellar Structure 
Program (see iFrost & Lattanziol Il996l and references therein 
for details) with AGB m ass loss based on the formula by 
IVassiliadis & WoodldI993h . The choice of mass loss is one of the 
most important uncertainties in the computation of stellar evolu- 
tion and it affects the nucleosynthesis by determining the stellar 
lifetime, which in turn changes the amount of TDU episodes as 
well as the timescale for the operation of HBB. Another impor- 
tant stellar model uncertainty is the treatment of convection and 
of convective borders. In our code we use the mixing length the- 
ory to parametrize convection with the free parameter a set to 
the value of L75. A different parametrization of convection and a 
different choice of the mass loss have been shown to have an im- 
portant im£actOTi_the_02eration of HBB and the resulting AGB 
yields (I Ventura & D'Anton a 2005a b). Furthermore, the treat- 
ment of convective boundaries affects the efficiency of the TDU 
(e.g. Frost & Lattanzio 1996). Here, we apply the Schwatzschild 
criterion and then extend the conve ctive zone t o the neutral bor- 
der in the same way as described in ' Lattanziol (Il986l) . 

The post-processing code calculates the abundances due to 
convective mixing and nuclear reaction rates for a large number 



of species. We use a nuclear network of 166 species from neu- 
trons and protons to S and from Fe to Nb. The missing species 
are accounted for by means of artificial neutron sinks. For one 
specific model, the 6 Mq Z = 0.02, we tested the results pro- 
duced using a full nuclear network of 320 species and no neutron 
sinks and found no differences (to within 0.1 dex) for [Rb/Fe], 
[Zr/Fe], and [ Rb/Zr]. The bulk of ou r 1285 reaction rates are 
from reaclib ('Thielemann et al.' 1 19861) . with updated tables as 
described by Lugaro et al. ( 2004 and iKarakas et al] (l2006h in 
particular for the '^^Ne{a, n)^^Mg r eaction. We furth er updated 
the neutron capture rates to those of iBao et"an (l2000l) . 

We have used here the same stellar models presented by 
iLugaro et al.l ( l2007h to compare to the composition of peculiar 
Stardust spinel grain OC2 and by Karakas et al. (2009) to com- 
pare to the observations of heavy elements in planetary nebulae 
of Type I. Th e Z ^ 0.008 and Z ^ 0.004 mo dels were also al- 
ready used bv lGarcfa-Hernandez et al.l (l2009l) to compare to the 
observations of intermediate-mass Rb-rich O-rich AGB stars in 
the Magellanic Clouds. 

A summary of the properties of the models is given in 
Table |2| including the mass (Mass) and metallicity (Z) of the 
model, the number of TPs (TPs), the maximum temperature 
achieved in the TPs {T'^^^), the maximum temperature achieved 
at the base of the convective envelope (T"^™"), the total mass 
dredged-up to the surface by the TDU (Mdred), the final enve- 
lope mass at the end of the computed evolution (Menv), if HBB 
is at work, and the final C/O and '^C/'^C ratios. 

The observations o f Garcfa-Hernandez et al.l (l2006h and 
iGarcia-Hernandez et al. I (l2007ir are of Galactic-disk AGB stars, 
and for this reason we concentrated on models of solar metal- 
licity Z = 0.02 for masses 3, 4, 5, 6, and 6.5 M^. The ini- 
tial abundances were taken from lAnders & Grevess 30989'). We 
also calculated a 6.5 Afo mode l using the solar abundances from 
lAsplundl dfOOSh and iLoddersI ([20031, which prescribe a solar 
metallicity of Z - 0.012. We further calculated low metallic- 
ity models of 6 M© Z = 0.008, 5 Mq with Z = 0.008, and 4 and 
5 Mq with Z = 0.004 which are suitable for comparison to the 
Magellanic^Cloud AGB stars. For the low-metallicity models we 
used the lAnders & Grevessd (Il989l) abundances scaled down to 
the appropriate Z. 

In four selected models we artificially included a PMZ in the 
top layers of the inters hell in the same way described in detail 
by Lugaro et aP (l2004l) . For the 3 Mq model we used a PMZ 
of size 2 X 10"^ Mq . For the 4, 5, and 6.5 models we used a 
PMZ of size 1 X 10^^ Mq . This smaller value is because the 
mass of the intershell is about one order of magnitude smaller 
in the higher than in the lower mass AGB stars (~ 10 Mq 
with respect to ~ 10"^ Mq), so we scaled the mass of the PMZ 
accordingly. The inclusion of a PMZ is the main uncertainty in 
models of the s process. We still do not know the mechanism 
and the features of the process driving the mixi ng, although some 
promising candidates have been proposed (see lBusso et al.l 19991 
for a discussion). Stellar models of intermediate-mass AGB stars 
indicate that any protons mixed into the He-intershell during the 
TDU would likely burn to '"^N before producing '''C owing to the 
high te mperat ures at the base of the envelope ( S iess et al. 20041: 
Her wTg||2004l) . This would inhibit the release of neutrons by the 
'^C neutron source. For this reason, the artificial inclusion of a 
PMZ into a model of an intermediate-mass AGB star should be 
treated with some caution. 
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Table 1. Ranges of the observed Rb and Zr abundances in AGB stars in the Galaxy and in the Magellanic Clouds. 





Putative low-mass AGB stars 


Putative intermediate-mass AGB stars 




Galaxy 


Magellanic Clouds" 


Galaxy 


Magellanic Clouds 


[Rb/Fe" ] 


-0.3 - -hO.9 


(-1.2- +0.4) 


+0.4 - +2.5 


+ 1.7 - +5.0 


[Zr^/Fe*] 


-0.3 - +1A 


(-0.3 - +0.70) 


0'' 


0'' 


[Rb/Zr] 


<0 


« 0) 


>0 


>0 


[Fe*/H] 


-1.3 - +0.1 


(-0.70 - -0.37'- ) 


set to 0' 


set to -0.3' and -0.7'' 




Lambert etal. (1995) 


Plez et al. (1993) 


Garcfa-Hernandez et al. (2006) 


Garcfa-Hernandez et al. (2009) 




Abia etal. (2001) 


Garcfa-Hernandez et al. (2009) 


Garcfa-Hernandez et al. (2007) 





Notes. As discussed in the text, these stars have high Li abundances and may represent a link between low-mass and in termediate-mas s AGB 
stars in the Magellanic Clouds. **' In some of the reported studies the average of Fe and Ni, M, is used instead of Fe. In lLambert et al.l ( Il995h 
Is/M] is reported as the average of several ^-process elements. We use this value as an indicator of the [Zr/Fe] ratio for stars from this study. 

The observational uncertainties allow for an upper limit of +0.5 and of +0.3 for Galactic and Magellanic Clouds intermediate-mass AGB stars, 
respectively. The metallicity for each star is not reported and the Rb and Zr abundances are derived using model atmospheres for the average 
abundances [M/H] set to 0, -0.3, and -0.7 in the Galaxy, Large and Small Magellanic Clouds, respectively. 

Table 2. Details of stellar models. 



Mass 


Z 


TPs" 


He 
(MK) 


2^ max 
bee 

(MK) 


(Mo) 


final Menv 
(Mq) 


HBB 


C/0 


"C/"C 


3 Mo 


0.02 


26 + 


302 


6.75 


0.081 


0.676 


No 


1.40 


117.47 


4Mq 


0.02 


18 + 2 


332 


22.7 


0.056 


0.958 


No 


0.99 


76.22 


5 Mo 


0.02 


24 + 5 


352 


64.5 


0.050 


1.500 


Yes 


0.77 


7.83 


6 Mo 


0.02 


438 + 5 


369 


83.1 


0.058 


1.791 


Yes 


0.38 


10.65 


6.5 Mq 


0.02 


40 + 7 


368 


86.5 


0.047 


1.507 


Yes 


0.34 


9.76 


6.5 Mq 


0.012 


52 + 2 


369 


90.0 


0.065 


1.389 


Yes 


0.75 


10.40 


5 Mo 


0.008 


59 + 4 


366 


80.8 


0.17 


1.795 


Yes 


0.93 


7.49 


6 Mo 


0.008 


69 + 5 


374 


89.6 


0.12 


1.197 


Yes 


1.41 


8.88 


4 Mo 


0.004 


30 + 2 


366 


66.3 


0.11 


1.333 


Yes 


3.18 


7.98 


5 Mo 


0.004 


83 + 3 


377 


84.4 


0.22 


1.141 


Yes 


2.59 


7.88 



Notes. +n refers to the number of synthetic TPs as described in Sect. l3.1l 



3.1. Late AGB-phase evolution 

Owing to convergence difficulties, the stellar evolution mod- 
els used as input into the post-processing calculations were not 
evolved through to the end of the superwind phased and into 
the post- AGB phase. That the stellar models have some remain- 
ing envelopes (see Table |2]i means that there could be further 
TPs and TDU events (Karakas & Lattanzio 2007). These latest 
stages of AGB star evolution are pivotal to the abundances un- 
der study here. Not only are the observed AGB stars in their late 
stages of evolution, but the freshly synthesized isotopes from 
the stellar interior are less diluted inside the rapidly decreasing 
envelope. To calculate the surface abundance evolution during 
the final stages of AGB evolution, we use a synthetic patch that 
we run after the nucleosynthesis code has finished. We estimate 
the remaining number of pulses based on the stellar parameters 
during the final time step that was calculated by the stellar evolu- 
tion code, as described in detail by Karakas & Lattanzio (2007). 
The number n of remaining pulses for each model are reported 
in Table [T] as +n in Column 3. For example, we estimate that 5 
and 7 more TP will occur for the 6 Mq and 6.5 Mq, Z = 0.02 
models, respectively. 

The most crucial component for our purposes in these cal- 
culations is the TDU efficiency as a function of the decreas- 
ing envelope mass. The TDU efficiency A is defined as the 
fraction of the mass that is dredged up over the mass by 
which the H-exhausted core had increased in the previous in- 



terpulse period, i.e., A = 



AMn 



The TDU efficiency in these 



' The superwind is the fin al stage of the AGB ph ase when very high 
mass loss, up to 10 '' Mo/yr dlben & Renziriilll983h may occur. 



final pha ses when the envelop e is decreasing could range from 



A ^ (i.S tranier oet allll997h to A 



0.9, at the last TP 



(IStancliffe & Jeffervl I2007L iKara kas 2010). In our case these 



(maximum) efficiencies are: A = 0.93, 0.95, and 0.94, for mod- 
els of 5, 6, and 6.5 Mq, respectively. This simple synthetic ex- 
tension of the evolution does not take into account the effect of 
HBB. However in all cases HBB has ceased in ffie detailed stel- 
lar evolution calculations. The Li abundance is not expected to 
vary significantly at this stage, owing to the depletion of the ^He 
fuel needed for Li production and the cessation of HBB. 



4. Results 

Figure |2] shows the maximum neutron density in the He in- 
tershell calculated for two of our models, a low-mass and a 
intermediate-mass AGB star, as a function of normalised model 
number (proxy for time). As anticipated, the maximum neu- 
tron density from the activation of the ^^Ne neutron source is 
up to two orders of magnitude higher in the intermediate-mass 
AGB than in the low-mass AGB model, where instead the '^C 
is the main neutron source. This results in positive and negative 
[Rb/Zr] in the two models, respectively. 

Our results for Rb, Zr, and Li are summarised in Table |3] 
and illustrated in Figs. |3] 111 E] and |6] As reference for future 
observations, in Table |4] we present the results for a set of key 
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Table 3. Final envelope abundances and ratios involving Rb, Zr, the isotopes of Rb, and Li for all our computed models. 



M 


Z PMZ mass (M©) 






[z^r/ rej 




|_ ISAj j JVUJ 




3Mo 


0.02 





26 


0.00 


0.00 


-0.01 


0.00 


1.51 


3 Mo 


0.02 


2 X 10-3 


26 


0.50 


1.00 


-0.49 


-0.11 


1.51 




0.02 





18 


0.01 


0.01 


0.00 


0.01 


1.30 


extended 


20 


0.03 


0.03 


0.00 


0.02 






0.02 


1 X 10" 


18 


0.57 


0.79 


-0.22 


0.08 


1.18 


extended 


20 


0.83 


1.21 


-0.38 


0.15 




5 Mo 


0.02 





24 


0.05 


0.01 


0.03 


0.07 


2.27 


extended 


29 


0.33 


0.16 


0.17 


0.37 




5 Mo 


0.02 


1 X 10^* 


24 


0.39 


0.40 


-0.01 


0.37 


2.37 


extended 


29 


0.80 


0.84 


-0.04 


0.63 




6 Mo 


0.02 





38 


0.21 


0.07 


0.14 


0.25 


1.06 


extended 


43 


0.73 


0.46 


0.27 


0.61 




6.5 Mo 


0.02 





40 


0.26 


0.09 


0.16 


0.30 


0.93 


extended 




1 C\A 


\j. 1 J 


n 'X 1 


n ^.Q 
u.oy 




6.5 Mo 


0.012 





52 


0.22 


0.07 


0.15 


0.25 


-0.49 


extended 


54 


0.69 


0.36 


0.33 


0.56 




6.5 Mo 


0.012 


1 X 10 


52 


0.62 


0.60 


0.02 


0.59 


-1.03 


extended 


54 


1.04 


1.04 


0.00 


0.76 




5 Mo 


0.008 





59 


0.71 


0.34 


0.37 


0.52 


-0.33 


extended 


63 


1.30 


1.02 


0.29 


0.76 




6 Mo 


0.008 





69 


0.90 


0.51 


0.39 


0.60 


-0.50 


extended 


74 


1.44 


1.10 


0.34 


0.77 




4 Mo 


0.004 





30 


0.80 


0.28 


0.52 


0.46 


2.98 


5 Mo 


0.004 





83 


1.12 


1.08 


0.04 


0.21 


-2.68 




extended 




86 


1.38 


1.45 


-0.08 


0.26 





Notes. Rows are alternately coloured for ease of reading. The rows labelled as extended report the final results obtained using the simple synthetic 
extension of the evolution described in Sect. 13. II 



Table 4. Same as Table[3]but for abundances and ratios involving C, N, O, Na, Mg and the isotopes of Mg for our computed models 
without inclusion of a PMZ. 



M 


Z 


[C/Fe] 


[N/Fe] 


[O/Fe] 


[Na/Fe] 


[Mg/Fe] 


[-■'Mg/-*Mg] 


pMg/^^Mg] 


3 Mo 


0.02 


0.49 


0.34 


-0.02 


0.19 


0.01 


0.03 


0.03 


4Mo 


0.02 


0.34 


0.38 


-0.02 


0.17 


0.03 


0.10 


0.10 


5 Mo 


0.02 


0.21 


0.57 


-0.04 


0.28 


0.04 


0.13 


0.18 


extended 


0.60 


0.56 


-0.05 


0.31 


0.14 


0.41 


0.50 


6 Mo 


0.02 


-0.14 


0.89 


-0.08 


0.85 


0.07 


0.22 


0.31 


extended 


0.46 


0.88 


-0.08 


0.85 


0.20 


0.51 


0.67 


6.5 Mq 


0.02 


-0.19 


0.87 


-0.09 


0.79 


0.06 


0.22 


0.30 


extended 


0.70 


0.86 


-0.09 


0.80 


0.32 


0.69 


0.89 


6.5 Mq 


0.012 


-0.16 


0.87 


-0.41 


0.64 


0.04 


0.46 


0.38 


5 Mo 


0.008 


0.19 


1.51 


-0.14 


1.52 


0.31 


0.59 


0.77 


6 Mo 


0.008 


0.26 


1.39 


-0.25 


1.04 


0.28 


0.79 


0.80 


4 Mo 


0.004 


0.88 


1.57 


-0.17 


0.93 


0.35 


0.68 


0.92 


5 Mo 


0.004 


0.56 


1.89 


-0.22 


1.78 


0.63 


0.83 


1.10 



Notes. Solar reference ratios are taken from lAnders & Grevessd ( Il989h . 



light elements (C, N, O, Na, and Mg) and the Mg isotopic ratios. 
Inclusion of a PMZ does not significantly affect the results for 
the light elements, except for the [Na/Fe] ratio in the 3 Mq and 
4 Mq models with Z = 0.02, which increases by roughly +0.14 
dex w ith the inclusion of the PMZ (see also Goriely & Mowlavj 
I2OOOI) . Same examples of the effect on the light elements of the 
synthetic extension of the models are also reported in Table |4] 
The yields of elements up S for a large range of stellar models, 
inclu ding those presented here, have been published in Karakas 
(120101 and re s ults fo r Zn, Ge, Se, Br, and Kr are published in 
iKarakaset aLl(l2009h . 

We predict an increase in [Rb/Fe] with increasing stellar 
mass, as is implied by the observations. This is because a higher 



stellar mass leads to (i) hotter temperatures during He-shell burn- 
ing and a stronger activation of the ^^Ne neutron source. This 
results in i-process enhancements of [Rb,Zr/Fe] in models more 
massive than 4 Mq even without the inclusion of the PMZ. Also, 
(ii) a higher number of TPs with high He-burning temperatures 
and followed by efficient TDU episodes (Table [TJ. However, 
intermediate-mass AGB models without synthetic extensions or 
PMZs produce [Rb/Fe] ratios that do not match the observed val- 
ues, while the predicted [Zr/Fe] ratios are within the observed 
range (< 0.5, Fig.©. 

We included a PMZ of 1 x lO""* Mq in our calculations for the 
5 Mq, Z = 0.02 and 6.5 Mq, Z = 0.012 models and found that 
[Rb/Fe] increases by a factor of 2 to 3. However, [Zr/Fe] also 
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normalised model number (proxy for time) 

Fig. 2. Maximum neutron density in the He intershell versus nor- 
malised model number (proxy for time) showing the last three 
TPs computed for a 3 Mg AGB star with a PMZ of 2 xlO"^^ 
Mq (red curve) and a 6.5 Mq model without a PMZ included 
(blue curve). The peaks are a result of the '^^Ne(a, n)^^Mg reac- 
tion activated during TPs. For the 3 Mq model the neutron flux 
during each TPs lasts for ~ 4 yr and the total time-integrated 
neutron flux (neutron exposure r) is ~ 0.02 mbarn"'. In the 6.5 
Mq model the neutron density reaches lO' cm"^ for ~ 3 yr and 
T ~ 0.2 mbarn '. The secondary bumps at ~ 10^ cm in the 
3 Mq model are a result of the ^^C{a, «)'^0 reaction, which is 
activated during the interpulse phase due to the PMZ. This acti- 
vati on lasts for ~ 30, 000 yr and results in r ~ 0.4 mbarn"' (see 
also lLugaro et al. 1^2003). The plateaux at ~ 10^ cm"^ represent 
neutrons released during radiative shell He burning in the inter- 
pulse periods. These neutron fluxes are insignificant compared 
to the others (with t < 0.02 mbarn"'). Furthermore, the mate- 
rial affected does not make it to the stellar surface as it remains 
buried in the C-O core. (Colour figure available online). 




[Zr/Fe] 



Fig. 3. Model results for the [Rb/Fe] against [Zr/Fe] for 
intermediate-mass AGB stars of solar metallicity without the in- 
clusion of a PMZ. The solid part of the lines are equivalent to a 
TDU efficiency of during the extended evolution. The dashed 
lines are the result of the synthetic extension when the TDU ef- 
ficiency is chosen to remain constant. The observed [Rb/Fe] and 
[Zr/Fe] ranges are indicated by the grey shaded area, note that 
[Rb/Fe] ranges up to 2.5 (+ 0.8, see Table [1]). (Colour figure 
available online). 




[Zr/Fe] 



Fig. 4. Results for selected intermediate-mass AGB models with 
and without the inclusion of PMZs. The solid lines represent 
model predictions without the inclusion of a PMZ, the dashed 
line models with a PMZ of 1 x 10""* Mq. The observed [Rb/Fe] 
and [Zr/Fe] ranges are indicated by the grey shaded area, note 
that [Rb/Fe] ranges up to 2.5 (+ 0.8, see Table[TJ. (Colour figure 
available online). 



increases by similar factors. This results in [Rb/Zr]~ 0, lower 
than found without the inclusion of a PMZ (Fig. HJi. A PMZ in 
intermediate-mass AGB stars may not provide a solution to the 
Zr and Rb abundances investigated here. 

In compariso n with the 5 Mq models presented by 
lAbia et all (1200 lb . we obtain lower production factors for Rb 
and Zr. This may be due to the different choice of mass-loss 
rate. lAbia et al. (2001) used the prescription by Reimers (1975|) 
with a choice of the free parameter y? - 10, while we used the 
prescription of'Vassilia dis & WoodI (Il993h . However, Table 1 of 
Straniero et al. (2000) indicates that the 5 Mq models used in 
lAbia et al.l (120011) experienced 23 TPs, very close to the number 
in our model. The difference must be due to a diff'erent '^C abun- 
dance in their PMZ. From comparison of the low-mass mod- 
els we derive that the results from our 3 Mq AGB model are 
in ag reement with the 1 .5 Mq model presented by lAbia et all 
(1200 1) when their standard choice of the amo unt of the ^^C neu- 
tron source is selected. For their 5 Mq model. IXbia et al. 
divide the amount of '^C by a factor of 10, while we divide the 
extent in mass of the PMZ by a factor of 20, which results in 
lower abundances. In any case, there is agreement in terms of the 
[Rb/Zr] ratio both specifically for the 5 Mq model, as well as in 
general, since our intermediate-mass AGB models all produce 
positive values ratios up to +0.4 (except in the 4 Mq Z - 0.02 
model with a PMZ included and the 5 Mq, Z - 0.004 model), 
similarly to the models of Abia_et al. ( 200 1 ) (except for one neg- 
ative value of -0.4 obtained for their specific standardx2 choice 
of the '^C amount). 

4.1. Synthetic Extension of the Models 

We extended our models to evolve through the superwind phase 
using the synthetic algorithm described in Sect. 13.11 The result- 
ing abundances are plotted in Fig. [3] (dashed parts of the lines). 
The solid parts of the curves in Fig. [3] are equivalent to assum- 
ing A = during the synthetic extension of the models, while 
the dashed parts of the curves show the results of extending 
the models using the same dredge up efficiency as during the 
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last computed pulse, i.e., the maximum possible dredge up effi- 
ciency. In this case, we find sufliciently high [Rb/Fe] at the final 
pulse for the 6 Mq and 6.5 Mq models to match the lowest range 
of the observations (within the observational error bar). We re- 
emphasise that this synthetic extension is a rough estimate. Even 
in this case we still cannot match the high end of the observed 
[Rb/Fe] values, i.e., > 1.8. These may be accounted for by AGB 
stars even more massive than 6.5 Mq. Another difficulty is that 
most likely HBB is shut off" in these late phases and hence it is 
not possible to keep the star oxygen rich as observed because 
carbon is carried to the surface via the TDU but it is not con- 
sumed by proton captures at the base of the envelope (Table|4]l. 

4.2. Reaction Rates 

Using the 6.5 Mq and 6 Mq Z - 0.02 models we performed sev- 
eral tests changing the neutron-capture cross sections involved 
in the production of Rb specifically with the aims of (i) inves- 
tigating the relative impact of the two branching points at ^^Kr 
and **^Rb on the production of Rb and (ii) finding the choice 
of neutron-capture rates, within their uncertainties, that would 
maximise the Rb production. The neutron-capture cross sections 
of the unstable isotopes '^^Kr and '^''Rb are evaluated only the- 
oretically, and the estimates vary roughly by a factor of two 
(see the www.kadonis.org database) thus we multiplied and di- 
vided them by a factor of 2. The neutron-capture cross sec- 
tion of **^Kr has been measured both via the activation tech- 
nique and the Time Of Flight (TOF) technique. The rate rec- 
ommended bv lBao et al.l(l2000.) is based on the activation mea- 
surements, however, the TOF measurements are typically higher 
(up to 70%) than the activation measurements. We made a test 
multiplying the ^^Kr(n, y)^^Kr cross section by a factor of 2. 
The neutron-capture cros s sections of **^Rb and '^^Rb have been 
recently re-measured by iHeil et all ([2008 ). who confirmed the 
values of the previous experiments with a high precision. We 
also made tests with some combinations of the above choices. 
The overall result is that the [Rb/Zr] ratio increased at most 
by 0.28 dex, when considering the synthetic extensions of the 
models and the following combination of test neutron-capture 
cross sections cr: o-[^^Kiin,jf^Ki]/2, cr[^''Kr(n,7)'*^Ki-]x2, and 
cr[**''Rb(n,7)**^Rb]x2. 

Given the surprising result that more Rb was produced when 
the efficiency of the ^^'Kr branching factor was smaller, we de- 
cided to test the differential impact of the **^Kr and the **^Rb 
branching points on the production of Rb using the 6 Mq 
Z = 0.02 model. Assuming instantaneous yS-decay of **^Kr we 
obtained [Rb/Fe] +0.14 higher because ['^^ Rb/Fe] decreased by 
-0.08 dex, but ['^''Rb/Fe] increased by -1-0.25 dex. We ascribe this 
result to the fact that when the ^^Kr branching point is open ^^Kr 
is produced, which has an extremely low Maxwellian-averaged 
neutron-capture cross section of 3.4 mbarn at a thermal energy 
of 30keV (Bee3[T99Tl) . and it accumulates at the expenses of 
^^Rb and ^^Rb. In fact, [^^Ki/Fe] decreased by -0.24 dex in this 
test case. Assuming instantaneous y6-decay of **^Rb we obtained 
instead [Rb/Fe] -0.08 lower as a consequence of a -0.25 dex 
lower [**^Rb/Fe]. 

Finally, again using the 6 Mq Z = 0.02 model, we checked 
the eff'ect of varying the ^^Ne(c!', n)^^Mg reaction rate. When 
varying the rate between the upper and lower limits described 
in|Karakas et al. (2006) we found variations in [Rb/Zr] of +0.16 
dex at m ost. When using the rate estimated in the NACRE com- 




z = 0.02 
z = 0.02 A 
= 0.008 
= 0.008 A ~ 
Z = 0.004 
Z = 0.004 X 
= 0.008 
= 0.008 A ~ 



1.5 



[Zr/Fe] 



Di lation (|Angulo et al. 1999i) . which is up 40% faster than that 
of lKarakas et al.T200( ) at typical He-burning temperatures, we 



found an increase of +0.21 dex in [Rb/Fe] and of +0.07 dex 



Fig. 5. Resulting [Rb/Fe] and [Zr/Fe] for low-metallicity 
intermediate-mass AGB models indicated by the labels. For 
comparison we included the 5 Mq Z = 0.02 model. As in Fig.|3] 
the solid part of the lines are equivalent to /I = during the ex- 
tended evolution. The dashed lines are the result of the synthetic 
extension when the TDU efficiency is chosen to remain constant. 
The observed [Rb/Fe] in five Rb-rich stars in the Magellanic 
Clouds ranges from +1.7 dex to +5.0 dex, outside of the range 
of the figure. (Colour figure available online). 



increase in [Zr/Fe], resulting in a +0.14 dex increase in the 
[Rb/Zr]. We have also computed selected models using the new 
-^Ne+Qf reaction rates from Iliadis et al. (2010). The main find- 
ing is that with the new rates both [Rb/Fe] and [Zr/Fe] can in- 
crease by up to roughly +0.2 dex, in the 6 Mq Z = 0.008 model, 
which is one of the models with the highest He-burning tem- 
peratures. For the other models the effect is smaller and basi- 
cally negligible at solar metallicity. The reasons are that (i) the 
new ^^Ne(Q', «)^^Mg rate is not very diff'erent from that used in 
our study and (ii) the fact that the revised ^^Ne(a, 'y)^^Mg rate is 
roughly an order of magnitude lower than that used here is not 
significant for the s process. This is because the amount of ^^Ne 
that burns is very small in any case: at most the ^^Ne intershell 
abundance decreases by a factor of two. Hence, that the two a- 
capture channels do not need to compete with each other for the 
22Ne. 



4.3. Low-metallicity models 

Fig. |5] shows that low-metallicity intermediate-mass AGB stars 
show significant enhancements of Rb and Zr, even without in- 
cluding the final synthetic pulses (see also Table O. This is 
marginally due to the temperature in TPs being slightly higher in 
these models, and mostly due to the effect of metallicity on the 
neutron-capture process itself. The neutron flux is proportional 
to ^^Ne/Z and the amount of ^-Ne in low-metallicity AGB stars 
has a strong primary component due to conversion of primary 
dredged-up '^C into ''^N, which is then converted into ^-Ne via 
double-Qf capture in the TPs. Hence the neutron flux increases 
with decreasing metallicity because the amount of material that 
can capture neutrons (both the light element poisons and the ^^Fe 
seeds) is smaller for lower Z. This results in a larger number of 
free neutrons and a stronger neutron flux. Also in these cases in- 
cluding TDU in the final pulses results in a significant increase 
in the surface [Rb/Fe], as well as of [Zr/Fe]. 
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Pulse 



Fig. 6. Li abundance, as log(Li/H) + 12, versus pulse number for 
our solar-metallicity models. (Colour figure available online). 



4.4. Lithium 

The evolution of Li at the stellar surface for our solar metallicity 
models is shown in Fig. |6] As the mass of the star increases so 
does the temperature at the base of the convective envelope dur- 
ing the AGB phase (Table |2]i. As the temperature increases the 
first effect on the Li abundance noticeable in the 4, 5, and 6 Mq 
models is the destruction of Li by HBB. This is because proton 
captures on ''Li are akeady efficient at temperatures around 20 
MK for the low density conditions of HBB. These temperatures 
are lower than those at which the "^He(^He,y)^Be reaction is ac- 
tivated (around 40 MK for HBB), which is responsible for the 
production of Li, after electron captures on ''Be. As the tempera- 
ture at the base of the envelope inreases as the star evolves along 
the AGB, the Cameron-Fowler mechanism for the Li production 
begins, leading in all the models to a Li-rich phase where the Li 
abundance reaches a maximum of log(Li/H) -H 12 ~ 4. The Li 
is mixed to hotter temperatures where it is eventually destroyed 
by HBB. In our models we find that the Li-rich phase lasts for 
X! 10^ years in a 6 Mq model. At some point the supply of ^He 
runs out and Li production by HBB ceases. At this point the star 
will no longer be observed as Li rich. Our results are in agree- 
men t with Li production via HBB presented by IVentura et al.l 
(1200 0) for models of metallicity 0.01. There, the Li abundance 
also first decreases, in the models of mass less than ~ 6 Mq, 
and then reaches a maximum close to log(Li/H) + 12 ~ 4 before 
declining to lower values. From Table |3]it can be noted that for 
a given initial mass, stars of lower metallicity have a lower final 
Li abundance. This is becasue they reach higher temperature at 
the base of the convective envelope (see Table |2]i and destroy Li 
more efficiently. 

In Fig.|2]we compare p redictions for selected models t o the 
observations of stars from Garcfa-Hernandez et al ] d2006h and 
iGarcia-Hernandez et al.l (l2007h for which both Rb and Li data 
are available. Shown are also the results for [Rb/Fe] for the syn- 
thetic extension of the models, which as noted above, do not 
include possible further decreases of Li. The observations show 
large variations in qualitative agreement with the large changes 
in Li seen in the models and due to the strong sensitivity of Li 
to the temperature of HBB. Th e stars in the Small Magellanic 
Clouds observed bv lPlez et al] (11993) have log(Li/H) + 12 be- 
tween H-1.9 and -1-3.5 and negative [Rb/Fe], which most likely 
reflects their initial composition. Our 5 MQmodel of Z = 0.02 
would provide a good fit to these observations because HBB is 




[Rb/Fe] 

Fig. 7. Li abundance versus [Rb/Fe] observed in intermediate- 
mass AGB stars (squares) and predicted by our models (solid 
lines: computed evolution, dashed lines: synthetic evolution as 
in Figs. [3] and nil. Different line colours indicate the following 
models: 5 M© Z = 0.02 with PMZ of size 0.0001 Mq (black); 
6 Mq Z = 0.02; (blue) 6.5 Mq Z ^ 0.02 (dark green); 6.5 Mq 
Z = 0.012 (red); 6.5 Mq Z ^ 0.012 with PMZ of 0.0001 Mq 
(light green). The error bars represent the maximum errors for 
each individual data point, though for sake of clarity they are 
included for only one data point. Arrows indicate upper limits. 
(Colour figure available online). 



producing Li and keeping the '^C/'^C ratio low, as observed, 
but without a significant increase in the [Rb/Fe] abundance. 
However, this metallicity is too high to be compared to the Small 
Magellanic Clouds. For models at the required Z = 0.004 metal- 
licity Rb production occurs already from the lowest masses at 
which HBB is activated (e.g., theA-MQ models at this metallicity 
in TableO. The stars observed by Plez et al. (1993) either (i) are 
all stars in the early phases of the thermally-pulsing AGB where 
Li production is active but the Rb abundance has yet to increase, 
or (ii) indicate that HBB needs to be more efficient and be acti- 
vated at masses where no substantial Rb production is expected. 
This may be achieved by a different description of convection, 
either via the Full Spectrum of Turbulen ce or equivalently by 
using a higher mixing length parameter a IVentura & D'Antona 
(2005a). 

5. Discussion and conclusions 

Our main conclusions can be summarised as follows: 

1. We predict an increase in the [Rb/Fe] ratio with increas- 
ing stellar mass. This is also suggested by the observations. 
However, our standard models are not able to quantitatively 
reproduce the high [Rb/Fe] values found in the observed 
intermediate-mass OH/IR stars. 

2. Inclusion of a partial mixing zone (PMZ) increases the final 
surface Rb content, but also increases the Zr content, result- 
ing in a [Zr/Fe] that is higher than the observed upper limit. 

3. The uncertainties in the neutron-capture cross sections rel- 
evant to Rb production do not introduce significant enough 
changes to our predicted [Rb/Fe] to match the observed data. 
Varying the ^^Ne(a;, n)^''Mg neutron source rate produces 
modest increases to the [Rb/Fe] ratio, but not enough to 
match the observed values. 
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4. In order to remove the whole envelope we synthetically ex- 
tend our models. By including these final pulses we find 
[Rb/Fe] ratios that are high enough to match the lower end of 
the observed data, taking into account the large error bars on 
the data. This relies on the TDU efficiency remaining as high 
during the final few TPs when the envelope mass is small as 
it was at the last computed TP. A better understanding of the 
TDU efficiency with receding envelope mass is essential and 
requires further investigation. However, one difficulty is that 
HBB is expected to be ineffective during this final stage of 
AGB evolution. This means the star will likely become car- 
bon rich if the TDU is efficient, while the observed stars are 
oxygen rich. 

5. More Rb (in absolute numbers) is produced as the metal- 
licity decreases. This reproduces the observed qualita- 
tive trend as derived from recent data for low-metallicity 
intermediate-mass AGB stars in the Magellanic Clouds 
dGarcia-Hernandez et al.l 1200 9*). However, also in this case 
the models cannot match ffie observations quantitatively. 

The main overall difficulty in matching the Rb and Zr ob- 
served abundances is that the s process cannot produce differ- 
ences larger than 0.5 dex between Zr and Rb as these two el- 
ements belong to the same ^-process peak. On the other hand, 
the observations show differences between these elements up to 
5 dex. Since the maximum [Rb/Zr] ratio produced by the s pro- 
cess is independent of the stellar models and their uncertainties, 
as it derives from the basic operation of the s process and the 
neutron-capture cross sections of the nuclei involved (which are 
know to the level of 5%), we need to turn to observational un- 
certainties to explain the mismatch for the [Rb/Zr] ratio. As dis- 
cussed in Sect. 12 and in iGarcia-Herna ndez et al. (2009) the so- 
lution to the problem may lie with an incomplete understanding 
of the atmospheres of luminous AGB stars. More realistic model 
atmospheres for intermediate-mass O-rich AGB stars (e.g., the 
inclusion of a circumstellar dust envelope and 3D hydrodynam- 
ical simulations) as well as NLTE calculations need to be devel- 
oped in order to shed some light on the observed discrepancy 
between the observations and the AGB nucleosynthesis theoret- 
ical predictions. 

Another possibility that needs to be carefully evaluated is 
whether the Rb and Zr abundances in these stars are affected 
by dust condensation. The condensation te mperature of Z r in a 
gas of solar-system composition is 1741 K (lLoddersll2003h . and 
it would be even higher in a gas where the abundance of Zr is 
enhanced by the s process. It is reasonable to expect that a sig- 
nificant fraction of the Zr abundance is removed from the gas by 
condensation into dust. On the other hand, the condensation tem- 
perature of Rb is much lower (800 K), and this element is not ex- 
pected to be removed from the gas. This effect could explain the 
observations, which samples the gas a nd not the dust, a round the 
star A similar idea was put forward bv lDinerstein et all (12006 ) to 
explain the non-detection of the ground-state fine-structure line 
of triply ionized Zr in several planetary nebulae known to have 
enhanced abundances of other light neutron-capture elements. 
To explain the most extreme observed [Rb/Zr] ratios requires 
that most of the Zr and none of the Rb has condensed into dust 
around ffiese stars. This possibility needs to be examined by de- 
tailed models of dust formation. 

If this "dust selection" effect is found to be significant, 
then stellar model uncertainties can also be invoked to play a 
role in the mismatch between data and predictions with regards 
to the [Rb/Fe] ratio. For example, the inclusion of a PMZ in 
intermediate-mass AGB stars may provide a solution to the Rb 



abundances investigated here, although as discussed in Sect. |3] 
this would be at odds with the current modelling of the PMZ 
in intermediate-mass AGB stars. Another complementary solu- 
tion may be found by increasing the number of TPs and TDUs 
episodes experienced by the stellar models. This number de- 
pends on the mass-loss rate of intermediate-m ass AGB stars, 
which is also highly uncertain. IVassiliadis & Wo od ( 1993) noted 
that there are optically bright long period variable stars with pe- 
riods of w 750 days that are probably intermediate-mass stars of 
~ 5 Mq. In order to prevent thei r mass-loss p rescription from re- 
moving these observed objects. IVassiliadis & Wood (1993) rec- 
ommend a modification to delay the onset of ffie superwind in 
stars of masses greater than 2.5 Mq. Note ffiat in the models 
presented here we have not used this modification, which would 
keep the mass-loss rate low for a longer time and would result 
in more TPs and TDU episodes. It would also keep the envelope 
mass relatively high, allowing for HBB to operate efficiently for 
more TPs. The outcome would likely be higher levels of Rb (and 
Zr) production in an O-rich envelope. Also, Rb production in 
AGB stars with masses greater than 6.5 Mq need to be com- 
puted. Finally, a higher temperature in the TP would lead to 
higher Rb production. This may be related to ffie penetration of 
flash-driven convection deeper into the core due to the inclusion 
of hydrodynamical overshoot. This would in crease the effi ciency 
of the ^-Ne neutron source as discussed bv iLugaro et alJ (l2003h . 
Detailed modelling is needed to explore the viability of these 
scenarios. 

Future observations of ffie isotopic composition of Mg in the 
stars considered here would provide further observational evi- 
dence that the -^Ne neutron source is activated in these stars as 
the neutron-rich isotopes of Mg, ^^Mg and ^^Mg are expected to 
be enhanced due to direct production via the ^^N e-Ha reactions 
(iKarakas & Lattanzidl2003t iKarakas et ai]l2006L and Table |4]i. 
Several more elements could be measured in AGB stars, and 
N would be particularly interesting to constrain HBB. Also the 
derivation of other i-process elements, e.g., Ba, La, Ce, and 
Pb, is possible and could help to further constrain the models. 
Though the predictions presented here involved a nuclear net- 
work up to Nb, we are currently working to extend these models 
up to Pb. AGB stars have complex atmospheres so which el- 
ements can actually be measured depends on the temperature, 
metallicity, and abundance for each individual star. 

Finally, the observations and models discussed here indicate 
that intermediate-mass AGB stars synthesize large amounts of 
Rb and that more Rb (in absolute numbers) is produced when 
lowering the metallicity. From this prediction, we extrapolate 
that intermediate-mass AGB model of metallicities typical of 
Globular Clusters (Z ~ 10 "^) will produce even larger abun- 
dances of Rb than presented here. If intermediate-mass AGB 
stars produced the heavy elements in Globular Cluster stars, then 
this would rule out these stars as candidates for ffie O, Na, Mg, 
and Al anomalies. This is because no correlation has been ob- 
served between the heavy elements and the light elements in 
Globular Cluster stars^ (Yong et al. 2006, 2008b a). A preHmi- 
nary study by iKarakas et al] ( f2010 ) appears to confirm our pre- 
dictions on the basis of two intermediate-mass AGB models, 
however, a larger set of detailed models at the relevant metallici- 
ties needs to be calculated. Furthermore, we also need to account 
for the data of Plez et al. ( 1993), which indicates that at metal- 
licities appropriate to the Magellanic Clouds there is a mass or 
evolutionary time when HBB is active but not the ^^Ne neutron 
source. A dedicated analysis of these stars will be the focus of a 
future work. 
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